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Abstract 

Walking technicolor theory attempts to realize electroweak symmetry breaking as the spon- 
taneous chiral symmetry breakdown caused by the gauge dynamics with slowly varying gauge 
coupling constant and large mass anomalous dimension. Many-flavor QCD is one of the candidates 
owning these features. We focus on the SU(3) gauge theory with ten flavors of massless fermions 
in the fundamental representation, and compute the gauge coupling constant in the Schrodinger 
functional scheme. Numerical simulation is performed with 0(a)-unimproved lattice action, and 
the continuum limit is taken in linear in lattice spacing. We observe evidence that this theory 
possesses an infrared fixed point. 
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I. INTRODUCTION 



While the standard model has been established through a number of experiments, un- 
natural hierarchies are present between the electroweak scale and the Planck scale and also 
among the fermion masses. Large Hadron Collider (LHC) is expected to give a new insight 
into these hierarchies. Among various new physics models proposed so far, Technicolor (TC) 
model [l| is one of the most attractive ones in these regards, as it does not require any fun- 
damental scalar particles, which cause the former hierarchy, and its extension, Extended TC 
model 3, has a possibility to generate the Yukawa hierarchy in a dynamical way. For recent 
review articles, see, for example, Refs. |3J. 

TC should be a strongly coupled vector-like gauge system, which triggers spontaneous 
chiral symmetry breaking (SxSB). It is widely known, however, that the simplest TC models 
obtained by rescaling ordinary QCD have already been ruled out by the S'-parameter {4] 
and FCNC jsj constraints. Refs. j^] suggested a series of TC models to circumvent the 
FCNC problem. Those TC models appeal to the gauge dynamics in which the effective 
gauge coupling constant runs slowly (i.e. "walks") at a relatively large value over a wide 
range of energy scale above the S^SB scale, and in which the chiral condensate gets large 
anomalous dimension. Such TC is called walking TC (WTC), and possible candidates have 
been enumerated through semi-quantitative analyses |7|. Since the dynamics that underlie 
WTC significantly differ from those of two or three-flavor QCD, the naive scaling argument in 
N c or Nf to estimate the S'-parameter would not work, and any quantitative predictions from 
WTC require solving nonperturbative dynamics explicitly. Lattice gauge theory provides a 
unique way to study this class of models from the first principles at present. 

Search for candidate theories of WTC is frequently linked to the iV/-dependent phase 
structure of the gauge theories. Let us take SU(3) gauge theory with Nf flavors of fermions in 
the fundamental representation as an example. According to the analysis of the perturbative 
/3-function, the system with large enough Nf (Nf > 16.5) is asymptotically non-free and 
trivial unless non-trivial ultraviolet fixed point exists. On the other hand, if Nf is sufficiently 
small (Nf < 3) the dynamics is QCD-like and thus in the chirally broken phase. It is believed 
that for the in-between Nf there exists a so-called conformal phase, where the coupling 
constant reaches an infrared fixed point (IRFP) without S^SB set in, but confinement may 
take place jsj. The range of Nf in which the conformal phase is realized is called conformal 
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window, and is represented by Nj nt < Nf < 16.5. It is then natural to speculate that the 
gauge dynamics slightly below Nj vlt exhibit the features required for WTC; slow running 
of the gauge coupling constant and S^SB. The first goal in the search for WTC is thus to 
identify Nf\ 

In the past years, many groups have used techniques of lattice simulations to search for 
Nj Tlt and/or WTC through hadron spectrum, eigenvalue distribution of Dirac operator, the 
behavior of running coupling constant, or renormalizati on g roup analysis of candidate theo- 
(9)]. For non-lattice studies, see, for example, Refs. 



ries 

many flavor QCD 



11] . Among various candidates, 
12Hl9|. sextet QCD 2CH25I]. and two-color adjoint QCD 26|-|3j| have been 



intensively studied. In this work, we focus on many flavor QCD with N c = 3 and fermions in 

running coupling constants were 
or eight- and twelve-flavor QCD using the Schrodinger functional (SF) scheme on 



calculated 
the lattice 



32J. They concluded that twelve-flavor QCD has an IRFP at g$ F ~ 5 while eight- 



flavor QCD does not. In practice, the study of the running coupling alone is supposed to be 
unable to fully exclude the possibility of a large IRFP because it requires lattice simulations 
at arbitrarily large coupling. Even worse, the unphysical, bulk first-order phase transition 



was found to occur in strong coupling regime of several gauge theories 



23 



341 ] . In such 

simulations, there exists an upper limit on the bare coupling at which lattice calculation is 



sensible. Nevertheless, because of the supports from the spectroscopy studies 



the conclusion in Ref. 12] that the eight-flavor QCD is QCD-like, i.e. Nf 1 * > 8, seems to 
be established nowadays. 



After the work of Ref. 12J , one group [15J has presented an evidence of the conformality 
of twelve-flavor QCD. The opposite conclusion, however, has also been reported by the 
other groups HQ. Therefore A^ nt < 12 is still under debate. Clearly the observed 
contradiction must be clarified before going further. While in the spectroscopy study of 
twelve-flavor QCD many sources of systematic uncertainties due to finite volume, taste 
breaking, chiral extrapolation, lack of continuum limit, etc., remain to be quantified, the 
calculation of the SF coupling constant of Ref. [12] appears, at present, to be less ambiguous. 
In such a circumstance, we are tempted to explore the dynamics of ten-flavor QCD. In this 
paper, we investigate, as a first step, the running coupling constant of ten-flavor QCD on 
the lattice to see whether it shows conformal behavior. We find that the running slows down 
and observe evidence that this theory possesses an infrared fixed point. 
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The paper is organized as follows. In sec. Ull we give remarks on how we identify IRFP 
on the lattice. Sec. II III summarizes the coefficients relevant to the perturbative calculation 
of the running coupling constant for later use. In sec. \IV\ the simulation setup including the 
definition of the running coupling constant in the Schrodinger functional scheme is presented. 
In sec. El we describe analysis method and present the numerical results. Sec. IVIl is devoted 
to the summary and outlook. 

II. REMARKS ON SEARCHING FOR IRFP ON THE LATTICE 

Since there exists a subtlety in proving the existence of IRFP with lattice gauge theory, in 
this section we briefly explain what is actually calculated and then give how to identify the 
existence of IRFP. Here we focus on the concept only. For further details of the calculational 
and analysis method that we take, see the following sections. 

In this work, we calculate the renormalized coupling constant in Schrodinger functional 
scheme at two different length scales, L and s ■ L. In practice, this is realized by repeating 
the calculation on two different volumes, Z 4 and (s-Z) 4 , at a common lattice bare coupling 
where I = L/a. We denote those couplings by u (or g 2 (g 2 , I)) and g 2 (g$, s ■ Z), respectively. 
Using those, we define the discrete beta function (DBF) by B(u, s, I) = 1/ g 2 {g^ s • I) — 1/u, 
where the rescaling factor s is arbitrary but is fixed to 2. If the DBF is free from lattice 
discretization errors, the sign of this quantity may directly tell whether the coupling constant 
increases or decreases against the scale change by s at the scale L, which is implicitly set 
by the value of u that we can choose. Since discretization errors do exist, however, we need 
to take the continuum limit. The a — > limit is taken for a fixed L, i.e. for a fixed u, by 
varying lattice spacing a. A series of the DBF thus obtained is then a function of /, and 
the I = L/a — > oo limit is expected to give the continuum limit. In summary, the DBF is 
constructed from a pair of lattice volumes (Z 4 , (s ■ Z) 4 ), and choice of larger Z results in the 
DBF closer to the continuum limit. 

In practice, lattice spacing is varied by changing the lattice bare coupling g%. If g 2 (gQ, h) 
turns out to be always larger than g 2 {g 2 , h) with Z2 > h, B(u, s, Z) < should hold for any 
Z and s > 1. In this case, the bare coupling at which g 2 (g^li) is equal to a fixed value u 
becomes small as lattice size l\ increases or one approaches the continuum limit. Thus the 
a — > limit is realized in the g 2 — > limit. This is the case for asymptotically free theories 
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with no IRFP such as ordinary QCD, and no subtlety is present. Even if an IRFP exists 
in such theories, the situation does not change as long as the input u is smaller than the 
IRFP, gf RFP . In other words, if the DBF extrapolated to I — > oo (or equivalently l/l — > 0) 
is negative, the limiting value is interpreted as the continuum limit and the possibility that 
an IRFP exists below u is excluded. 

When the DBF extrapolated to I — > oo is positive, interpretation of numerical results 
becomes ambiguous. In this case, in the vicinity of l/l = 0, g 2 (gLs ■ I) < g 2 (gQ,l), i.e. 
B(u,s,l) > 0. Indeed, it happens below = 4.4 in Fig. 4 of Ref. 13|, for example. Then, 
one may expect that the I — > oo limit is realized by g 2 — > oo on first sight. However, 
recalling 4 theory, this expectation turns out to be too naive. In 4 theory, the continuum 
limit exists only in the trivial case unless the theory possesses a non-trivial UV fixed point. 
Since the situation is similar to this case, the most plausible interpretation is that, when 
u > g 2 RFP , the continuum limit does not exist unless a nontrivial UV fixed point exists. 
Since no nontrivial UV fixed points has been established so far, it is not suitable to call the 
extrapolated value the continuum limit when it is positive. Nevertheless, we can still infer 
that u > gjjtFP because no other possibility remains. 

We investigate the sign of the DBF, starting with the weak coupling regime u ~ 1 where 
the perturbative calculation is reliable and predicts a negative value. We keep monitoring 
the sign of the DBF with increasing u. The identification of the IRFP is then made by sign- 
flip of the DBF extrapolated to / — > oo. Notice that, when the extrapolated value is positive, 
the extrapolation does not make sense and hence we do not insist that the continuum limit 
is determined. 



III. PERTURBATIVE ANALYSIS 

We start with defining the /3 function of an effective gauge coupling constant in a mass- 
independent renormalization scheme, which should have the following expansion in the per- 
turbative regime 

P(g 2 (L)) = L = b, g\L) + b 2 g\L) + b, g\L) + 6 4 g 10 (L) + ■ • • , (1) 



5 



where L denotes a length scale. The first two coefficients on the right hand side are scheme- 
independent, and given by 



(4tt) 
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-N 



(4tt) 



102 - —N f 

3 i 



(2) 



The remaining coefficients are scheme-dependent and known only in the limited schemes and 
orders. The third coefficient takes the following form in the Schrodinger functional scheme; 
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2tt 
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where b^ s is a coefficient in the MS scheme, 
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(3) 



(4) 



and the calculable quantities c 2 and Cg depend on the spatial boundary condition imposed 
on the fermion fields in the SF setup, i.e so-called 9. Those for 9 = n/5 and c 2 for 9 = are 
known to be 35 ] 



6>=tt/5 
C 2 

0=tt/5 



J=0 



1.25563 + 0.039863 x iV 



0=7r/5v2 



1.197(10) + 0.140(6) x N f - 0.0330(2) x N 



2 

f 



1.25563 + 0.022504 x N 



(5) 
(6) 
(7) 



but C3 =0 has not been calculated yet. Although 9 = is chosen in our simulation as described 
in sec. IIVt the coefficients for 9 = ir/5 are used only to see the situation of conformal windows 
inferred just from the perturbative analysis, and the potential size of difference between the 
two- and three-loop calculations. 

The perturbative estimates of the infrared fixed point (IRFP) for SU(3) gauge theory 
with Nf flavors of fundamental fermion are summarized in Tab. HI We note that in the 
three-loop perturbative analysis the existence of IRFP is determined only by the sign of 63, 
which is always negative for the range of Nf shown in Tab. [H Therefore, the existence of 
IRFP as well as its value may be unstable against including higher orders. Nevertheless, for 
Nf > 14 the difference between the two- and three-loop results is reasonably small, and one 
may expect that higher order corrections do not spoil the existence of IRFP or even do not 
change its value by much for such a large Nf. 

According to the analysis based on Schwinger- Dyson equation, S^SB is expected to occur 
when the running coupling constant reaches g 2 ~ it 2 in SU(3) gauge theories [36j. In spite 
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Nf 


4 


6 


8 


10 


12 


14 


16 


two-loop universal 










27.74 


9.47 


3.49 


0.52 


three-loop SF with 9 = 


vr/5 


43.36 


23.75 


15.52 


9.45 


5.18 


2.43 


0.47 



TABLE I: The perturbative IRFP obtained from the two-loop universal and the three-loop SF 
scheme analyses. 

of the scheme-dependence of the running coupling constant and the value of IRFP, those 
results motivate us to speculate that ten-flavor QCD may exhibit strongly coupled walking 
dynamics, and thus deserves full nonperturbative calculation. 



IV. SIMULATION DETAILS 
A. Schrodinger functional 



We employ the Schrodinger functional (SF) method |32| to study the scale dependence 
of the running coupling constant. Unimproved Wilson fermion action and the standard 
plaquette gauge action are used without any boundary counter terms as described below. 

The SF on the lattice is defined on a four dimensional hypercubic lattice with a volume 
(L/a) 3 x (T/a) in the cylindrical geometry. Throughout this work, the temporal extent T/a 
is chosen to be equal to the spatial one L/a. Periodic boundary condition in the spatial 
directions with vanishing phase factor (9 = 0) and Dirichlet one in the temporal direction 
are imposed for both gauge (U^(x)) and fermion (ip(x) and ip(x)) fields. The boundary 
values for gauge and fermion fields are represented by three-by-three color matrices, C and 
C, and spinors, p, p', p and p', respectively. The partition function of this system is given 
by 

Z SF (C, p', p' ; C, p, p) = e - r «W ^ p) = J D[U, tf>, ^ e -m^,c,c', P , P ',p,p>] : ( 8 ) 
where T is the effective action, and 

S[U, ^ C, C, p, p', p, p'] = S g [U, C, C] + S q [U, 1>, A P, P\ P, A (9) 
For the pure gauge part, we employ the plaquette action, 

3 3 

S g [U, C,C'] = ^J2J2Y1 ^uw^xo) Tr [1 - P^(x)} , (10) 

x fi=0 i/=0 



where j3 = denotes the inverse of the bare coupling constant, 5^=0 when fi = v 

otherwise 1, and P fJ-jU (x) denotes a lxl Wilson loop on the \i-v plane starting and ending 
at x. The spatial link variables on the boundaries, the hypersurfaces at x — and L/a, are 
all set to the diagonal, constant SU(3) matrices as 



/ 



U k (x)\ xo =o 



exp [C] , C 



to 



U k (x) 



\xo=L/a 



exp [C'\ 



n> - m 



77-f ^ 

-\r] 
-h + fy 

77 — 7T 

h + f 



(11) 



V 



1 - 1 27r 



(12) 



where k — 1, 2, 3, and 77 is parameterizing the gauge boundary fields. The weight w^ u (x ) 
in eq. f lTUl) is given by 



c f for (t = or t — (L/a) — 1) and (// or i/=0) 

for (i = (L/a)) and (/x or z/=0) 

for (t — or t = (L/a)) and (/i 7^0 and 1/ 7^0) 

1 for all the other cases 



(13) 



By tuning q, 0(a) errors induced from the boundaries in the time direction can be removed 
perturbatively, but in this work we simply take its tree level values, c t = 1. With this setup, 
the value of c s can be arbitrarily chosen because the spatial plaquettes on the boundaries 
do not contribute to the action. We thus set c s = 0. 

The fermion fields are described by the unimproved Wilson fermion action, 

S q [U,il>,$\ = N f J2$(x)D(x,y;U)i>(y) = N f ^^(z)^^; U)^\y), (14) 

x,y x,y 

D Xat (x,y;U) = 5 xy - kJ^ {(1 - Ta») u ^) 5 *+fi,v + i 1 + Tm) u l( x ~ fi s *-fi,v} » ( i5 ) 



where 



^ lat (x) 



1 



ip(x), ip lat (x 



1 



1>(x), D lat (x,y;U) = 2 K D(x,y;U). (16) 



V2k V2k 

The hopping parameter k is related to the bare mass tuq through 2 k = l/(arrio + 4). The 
dynamical degrees of freedom of the fermion field ip(x) and anti-fermion fields ip(x) reside 
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on the lattice sites x with < xq < T. On both boundaries (xq = and T), the half 
of the Dirac components are set to zero and the remaining components are fixed to some 
prescribed values, p, p, p' and p', as 



P + ip{x)\ Xo=0 = p(x), 


P-ip{x) 


L =o 


= o, 


(17) 


P-4j(x)\ Xo=T = p'(x), 


P + ^{x) 


\x =T 


= o, 


(18) 


^(x)P-| xo=0 = p(x), 


4j(x)P + 


Lo=0 


= 0, 


(19) 


$(x)P + \ xo=T = jf(x), 


4>(x)P^ 


L =T 


= 0, 


(20) 



where P± = (1 ± 7o)/2. In this work, the boundary values for the fermion fields are set to 



p = p' = p = p' = 0. (21) 



B. Definition of the running coupling 



With the gauge boundary conditions ffTTj) and ffT2l . the absolute minimum of the action 
is given by a color-electric background field denoted by B(x). Then, the effective action can 
be defined clS 8b function of B by 



T[B] = -\nZ SF (C\p',p';C,p,p), 

which has the following perturbative expansion in the bare coupling constant, 

1 



r 



9o 



T + T 1 + 0(g t 



o 2 ), 



and, in particular, the lowest-order term 

T = [g 2 S 9 [B]} 



90=0 



(22) 



(23) 



(24) 



is exactly the classical action of the induced background field. The SF scheme coupling is 
then defined in the massless limit for fermions by 



dT 



k 



r)=0 



aipial iy 



(25) 



v=0 9h(glJ = L/a) dr] 
where the normalization constant k is determined such that g§ F = g\ holds in the leading 
order of the perturbative expansion, and is found to be 



k 



9To 
dr} 



rj=0 



12 Q [sin (2 7 ) + sin (7)] = k with 7 = | 



(26) 



Because of the absence of the clover term, only the ^-derivative of the gauge action con- 
tributes to 1/<7sfG7o) 0- 

C. Parameters 

The simulation was performed on the lattice sizes of / 4 = (L/a) A = 4 4 , 6 4 , 8 4 , 12 4 , and 
16 4 with the inverse of bare gauge coupling constant (3 = Q/g^ in the range, 4.4 < < 96.0. 
However, the data from / = 4 lattices are not used in the following analysis because it was 
found that they have large discretization errors. We calculated the SF coupling on 18 4 lattice 
with a single (3 (/3=4.55), and the result is used to check the scaling violation at a specific 
value of <7§ F . 

The algorithm to generate the gauge configuration follows the standard HMC with five 
pseudo-fermion fields introduced to simulate the ten flavors of dynamical fermions. The 
numerical simulations were carried out on several different architectures including GPGPU, 
PC cluster and supercomputers. In order to achieve high performance on each architecture, 
the HMC code, especially the fermion solver part, were optimized depending on each ar- 
chitecture. In particular, mixed precision solver using multiple GPUs enables us to obtain 
high statistics on g$ F at I 4 = 12 4 and 16 4 [37|. Acceptance ratio is kept to around 80 % by 
adjusting the molecular dynamics step size (St). 

Since the Wilson fermion explicitly breaks chiral symmetry, the value of k is tuned, for 
every pair of (/3, L/a), to its critical value k c realizing the massless fermion by monitoring 
the corresponding PCAC mass. The values of /3, k, the number of trajectories, St and the 
results for I = L/a=6, 8, 12, 16, and 18 lattices are tabulated in Tabs. iHWIt respectively. 

D. Comment on O(a)-unimprovement 

In our pilot study, we employed the 0(a)-improved fermion action with the perturbatively 
determined counter terms. With this setup, we encountered a sudden change of the plaquette 
and the PCAC mass at 1=6 and /3=3.6 when k was decreased from 0.1517, and we could not 
realize the vanishing PCAC mass. The expected SF coupling constant is about 3 ~ 4 there. 



The same phenomenon also occurs on / = 4 lattices at almost the same value o: 
constant. Since the observed behavior looks similar to those reported in Refs. 



bare coupling 



23 



33 



34j, we 
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infer that this is a bulk, first order phase transition. In order to cover the region g$ F ~ O(10), 
we omitted any 0(a) improvements. Thus the leading discretization error in our result is 
linear in lattice spacing. 

Even without 0(a) improvements, the bulk, first order phase transition is observed for 
(3 = 6/<7q ~ 4.4. However, this time it happens at the renormalized coupling constant greater 
than the 0(a)-improved case, typically g| F ~ 0(10). Since this bulk phase transition is 
considered as a lattice artifact, whenever this happens we discard the gauge configurations 
at such /3. Thus the position of the critical /3 (~ 4.4) sets the lower limit on our exploration 
of/?. 



V. ANALYSIS METHOD AND RESULTS 



A. Raw data 

The SF coupling constant (g$ F ) and the PCAC mass (M) obtained on each ((3, k, I) are 
shown in Tabs. iHllVIl g^/gsp is plotted as a function of the bare coupling constant g^ in 
Fig. [TJ The figure shows that g% F increases with I — L/a at a fixed g 2 , but the change 
between the data from / = 12 and Z = 16 is tiny. For later use, we fit the data of g$/ g$ F 
to an interpolating formula as a function of the bare coupling constant g^. Among various 
functional forms we examined, the following form 

9l = 1 ~ ai,i 9o / 2? x 

9h(9l " 1 + Pi,i x g\ + En= 2 ai,n x g 2 n ' 

turned out to give the minimum % 2 /dof for a fixed number of free parameters, N . We thus 
employ eq. (ETj) . In eq. (|2"T|) . pxj is the /-dependent coefficient and we have calculated them 
perturbatively in the SF scheme 



Pi,i 



0.4477107831 


for I = 


6 


0.4624813408 


for / = 


8 


0.4756888260 


for / = 


12 


0.4833079203 


for / = 


16 


0.4864767958 


for / = 


18 



(28) 



The other coefficients a; jn 's are determined for each / independently. We optimize the degree 
of polynomial N in the denominator of eq. (I27p by monitoring x 2 /dof, and take N = 5 for 

11 



I = 6 and 12, and N = 4 for I = 8 and 16. Tab. IVHI shows the fit results for the coefficients 
in eq. f )27|) . The fit results are also shown as the region sandwiched by a pair of solid curves 
in Fig. [TJ 
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FIG. 1: g 2 . dependence of So/fsF f° r ' = L/a=6, 8, 12 and 16. The right panel magnifies the region 
offlg€[l.l, 1.40]. 

Hereafter we denote the SF coupling obtained at a bare coupling constant g 2 and at a 
lattice length of I by <7sf(<7ojO an d its continuum counterpart by g§ F (L). 



B. Discrete /? function 

In order to see the scale dependence of the SF coupling constant, we analyze the discrete 



23| . The whole procedure is described below. 



/3 function (DBF) introduced in Refs. [20 

First, we choose an initial value of the running coupling constant, denoted by u. This 
implicitly sets the initial length scale Lq through gg F (Lo) = u. Using the interpolating 
formula (1271) for the lattice size I (= L/a), the bare coupling constant g$ is numerically 
obtained by solving the equation S^p^o 2 ' = u - I is identified with L /a, so that the 
lattice spacing at g$ 2 is found to be a(gQ 2 ,l) = Lq/1. Now we choose a rescaling factor, s. 
The lattice step scaling function T,q(u, s, /) is then defined as the SF coupling for I' = s ■ I at 
the same bare coupling #q 2 , i.e. 

S o (n,s,0 = 9lM 2 ,s-l)\ g2sAgfl)=u . (29) 

The meaning of the subscript "0" becomes clear soon. Of course, both I and s ■ I must be 
equal to one of 6, 8, 12 and 16, and hence the possible values for the rescaling factor s are 

12 



P_ 

96.0000 
96.0000 



K 



Trajs. 



plq. 5t 



Acc. 



9sf 



0.1267030 39,700 0.979268 
0.1267070 49,900 0.979267 



0.000002 
0.000002 



0.0076 0.827 
0.0076 0.826 



0.002 
0.002 



0.06431 
0.06428 



0.00006 
0.00005 



M 

0.00012(0.00003) 
-0.00004(0.00003) 



48.0000 0.1276060 40,100 0.958852 
48.0000 0.1276100 41,100 0.958846 



0.000005 
0.000004 



0.0098 0.857 
0.0098 0.857 



0.002 
0.002 



0.13221 
0.13209 



0.00010 
0.00010 



-0.00013(0.00002) 
-0.00016(0.00002) 



24.0000 0.1295180 24,700 0.917566 
24.0000 0.1295200 60,300 0.917562 



0.000009 
0.000005 



0.0149 0.848 
0.0152 0.838 



0.002 
0.002 



0.28079 
0.28086 



0.00015 
0.00010 



0.00006(0.00002) 
0.00006(0.00001) 



12.0000 0.1339640 48,700 0.833056 



0.000014 



0.0250 0.812 



0.002 



0.64450 



0.00076 



-0.00004(0.00005) 



9.6000 0.1365680 160,300 0.789765 



0.000012 



0.0256 0.826 



0.001 



0.87189 



0.00068 



0.00002(0.00003) 



7.4000 0.1410690 120,500 0.724148 



0.000015 



0.0270 0.854 



0.001 



1.30413 



0.00194 



0.00006(0.00006) 



6.8000 0.1430520 120,300 0.698517 



0.000018 



0.0270 0.870 



0.001 



1.51024 



0.00244 



-0.00025(0.00008) 



6.3000 0.1451400 17,400 0.673231 



0.000076 



0.0333 0.817 



0.001 



1.74788 



0.00505 



0.00015(0.00019) 



6.0000 0.1466380 33,600 0.655993 
6.0000 0.1466410 80,300 0.655981 



0.000034 
0.000029 



0.0333 0.837 
0.0333 0.833 



0.002 
0.001 



1.93684 
1.93605 



0.00691 
0.00362 



0.00044(0.00021) 
0.00004(0.00010) 



5.5000 0.1497590 50,300 0.622923 
5.5000 0.1497610 36,000 0.622942 
5.5000 0.1497620 140,300 0.622977 



0.000025 
0.000027 
0.000023 



0.0370 0.817 
0.0357 0.827 
0.0357 0.831 



0.002 
0.002 
0.001 



2.38340 
2.36232 
2.37542 



0.01092 
0.00634 
0.00963 



0.00042(0.00020) 
-0.00018(0.00022) 
0.00023(0.00014) 



5.2000 0.1521330 220,300 0.600097 



0.000019 



0.0380 0.812 



0.001 



2.80668 



0.01246 



-0.00015(0.00014) 



5.0000 0.1539800 59,900 0.583463 



0.000049 



0.0400 0.806 



0.002 



3.28837 



0.06618 



-0.00005(0.00046) 



4.6000 0.1585140 33,800 0.545776 
4.6000 0.1585150 150,000 0.545680 



0.000055 
0.000041 



0.0400 0.813 
0.0400 0.813 



0.002 
0.001 



5.47008 
5.41263 



0.13064 
0.09891 



0.00092(0.00043) 
0.00123(0.00042) 



4.5000 0.1599020 100,300 0.535280 
4.5000 0.1599030 100,300 0.535305 



0.000061 
0.000066 



0.0400 0.813 
0.0400 0.813 



0.001 
0.002 



7.02516 
6.70575 



0.24479 
0.19622 



0.00111(0.00069) 
0.00033(0.00061) 



4.4215 0.1610680 105,900 0.526537 
4.4215 0.1610820 92,400 0.526692 



0.000087 
0.000066 



0.0385 0.825 
0.0385 0.826 



0.001 
0.001 



8.90139 



0.36944 
0.32355 



0.00238(0.00097) 
0.00073(0.00075) 



4.4000 0.1614210 249,500 0.524331 
4.4000 0.1614220 182,500 0.524342 
4.4000 0.1614230 250,500 0.524387 



0.000060 
0.000091 
0.000062 



0.0400 0.811 
0.0400 0.812 
0.0400 0.811 



0.001 
0.001 
0.001 



9.60163 



0.19661 



10.17980(0.33990) 
10.07713(0.25379) 



0.00051(0.00050) 
0.00119(0.00073) 
0.00049(0.00053) 



TABLE II: Simulation parameters and results obtained at L/a=6. 



13 



(3 k Trajs. plq. 5t Acc. <7g F M 



96.0000 0.1263270 22,500 0.979420(0.000002) 0.0056 0.811(0.004) 0.06434(0.00004) 0.00001(0.00001) 

48.0000 0.1272250 18,300 0.958843(0.000003) 0.0100 0.818(0.007) 0.13247(0.00017) 0.00002(0.00002) 

24.0000 0.1291450 42,300 0.917260(0.000004) 0.0125 0.804(0.003) 0.28282(0.00023) -0.00004(0.00002) 

12.0000 0.1335850 68,500 0.832266(0.000007) 0.0167 0.828(0.005) 0.65380(0.00071) -0.00010(0.00003) 

9.6000 0.1361800 21,820 0.788830(0.000013) 0.0200 0.828(0.006) 0.88751(0.00287) -0.00008(0.00005) 

7.4000 0.1406600 63,330 0.723081(0.000010) 0.0250 0.818(0.003) 1.34182(0.00417) -0.00004(0.00016) 

6.8000 0.1426200 41,500 0.697409(0.000013) 0.0250 0.797(0.002) 1.56232(0.00662) 0.00012(0.00011) 

6.3000 0.1447000 28,000 0.672208(0.000021) 0.0250 0.816(0.003) 1.81987(0.01036) -0.00034(0.00014) 

6.0000 0.1462000 47,000 0.654999(0.000012) 0.0250 0.820(0.003) 2.01248(0.01258) -0.00042(0.00011) 

5.5000 0.1492700 35,900 0.622016(0.000021) 0.0286 0.797(0.003) 2.48139(0.01969) -0.00021(0.00015) 

5.0000 0.1533600 27,900 0.582458(0.000038) 0.0250 0.825(0.004) 3.46930(0.07238) 0.00094(0.00034) 

4.8000 0.1554270 114,500 0.564464(0.000020) 0.0250 0.860(0.001) 4.35348(0.09845) 0.00026(0.00024) 

4.7000 0.1565500 35,400 0.554789(0.000040) 0.0256 0.854(0.002) 4.87595(0.21035) 0.00027(0.00051) 

4.6200 0.1575500 86,300 0.546856(0.000030) 0.0312 0.783(0.001) 6.23744(0.25321) -0.00023(0.00033) 

4.6000 0.1577800 149,300 0.544695(0.000027) 0.0250 0.852(0.002) 6.01108(0.17093) -0.00007(0.00028) 

4.5500 0.1584200 24,500 0.539428(0.000090) 0.0278 0.833(0.003) 6.92022(0.46491) 0.00087(0.00088) 

4.5500 0.1584270 93,300 0.539336(0.000033) 0.0278 0.831(0.002) 6.99432(0.31873) 0.00135(0.00041) 

4.5500 0.1584500 25,700 0.539683(0.000064) 0.0278 0.832(0.004) 6.74187(0.46970) -0.00163(0.00071) 

4.5200 0.1588500 56,570 0.536316(0.000058) 0.0278 0.827(0.002) 8.28029(0.57687) -0.00010(0.00059) 

4.5000 0.1591300 107,100 0.534108(0.000036) 0.0250 0.859(0.001) 8.40630(0.37369) -0.00007(0.00038) 

4.4800 0.1594000 41,555 0.531781(0.000085) 0.0250 0.827(0.002) 8.57214(0.57202) 0.00027(0.00070) 

4.4215 0.1602640 160,900 0.525143(0.000050) 0.0263 0.837(0.001) 12.21877(0.49625) -0.00012(0.00041) 

4.4215 0.1602700 127,500 0.525149(0.000058) 0.0250 0.861(0.001) 12.62365(0.68980) -0.00059(0.00048) 

4.4200 0.1602700 29,700 0.524651(0.000132) 0.0278 0.828(0.002) 13.15085(0.99774) 0.00214(0.00075) 

4.4000 0.1606000 229,500 0.522502(0.000057) 0.0278 0.819(0.002) 15.00764(0.69115) 0.00020(0.00042) 

TABLE III: Simulation parameters and results obtained at L/a=8. 
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U000012) 


0.0182 0.818(0.002) 


7.29885(0.38589) 


0.00075(0.00017) 


4.5500 0.1578500 


57,260 


0.539025 


;o.oooo2i) 


0.0192 0.802(0.002) 


10.15231(1.11827) 


0.00125(0.00027) 


4.5000 0.1585500 104,570 0.534014 


;0.000049) 


0.0250 0.701(0.007) 


13.03915(1.33994) 


-0.00237(0.00039) 



TABLE IV: Simulation parameters and results obtained at L/a=12. 

limited. The difference between S (u, s,l) and u gives the scale dependence through the 
scale change from L to s ■ L, up to lattice artifacts. 

Since the raw data of 1/c7sf(s , 5 fluctuate around zero in the strong coupling region, 
converting from 1 /^sfCs'o' ^° 9sf(9o^) sometimes induces huge statistical uncertainty. To 



avoid this we treat the inverse coupling constant, l/g^ F (gQ,l), directly. Then, 
scale dependence of the inverse coupling constant, we introduce the lattice DBF 

1 1 



to see the 



20 



B (u,s,l) 



23|by 
(30) 



S (m, s,/) u 

We calculate the continuum limit of this function for various initial values of the coupling 
constant, u. If the sign of the DBF in the continuum limit turns out to flip at a certain 
renormalized coupling constant u, it indicates the existence of IRFP. 
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(3 k Trajs. plq. 5t Acc. <7g F M 



24.0000 0.1288000 22,890 0.916505(0.000002) 


0.0067 0.796(0.003) 


0.28368(0.00057) 


0.00011(0.00002) 


12.0000 0.1332590 40,950 0.830799(0.000003) 


0.0091 0.718(0.011) 


0.66384(0.00293) 


-0.00016(0.00002) 


9.6000 0.1358600 30,900 0.787096(0.000002) 


0.0080 0.807(0.003) 


0.90538(0.00570) 


-0.00023(0.00003) 


7.4000 0.1403250 62,300 0.721190(0.000003) 


0.0111 0.812(0.002) 


1.39094(0.00834) 


-0.00002(0.00003) 


6.8000 0.1422900 39,796 0.695613(0.000004) 


0.0133 0.787(0.002) 


1.63562(0.01649) 


-0.00042(0.00005) 


6.3000 0.1443400 69,000 0.670503(0.000004) 


0.0133 0.798(0.002) 


1.91412(0.01628) 


-0.00036(0.00004) 


6.0000 0.1457950 18,900 0.653363(0.000007) 


0.0156 0.712(0.004) 


2.12147(0.03887) 


0.00043(0.00010) 


5.5000 0.1488500 50,330 0.620911(0.000007) 


0.0143 0.782(0.002) 


2.67936(0.03897) 


-0.00047(0.00011) 


5.0800 0.1521310 23,760 0.588803(0.000007) 


0.0139 0.804(0.003) 


3.24742(0.07271) 


-0.00003(0.00012) 


5.0000 0.1528550 71,954 0.582121(0.000004) 


0.0143 0.797(0.002) 


3.86709(0.12622) 


-0.00004(0.00009) 


4.8000 0.1548310 46,000 0.564445(0.000008) 


0.0156 0.755(0.002) 


5.72911(0.49013) 


0.00003(0.00014) 


4.6000 0.1570500 83,705 0.544764(0.000008) 


0.0143 0.794(0.001) 


8.21243(0.63114) 


0.00128(0.00012) 


4.5500 0.1576750 107,069 0.539609(0.000010) 


0.0139 0.809(0.002) 


10.81452(0.80073) 


0.00011(0.00011) 


4.5200 0.1580650 42,400 0.536387(0.000021) 


0.0156 0.754(0.002) 


17.34193(3.72829) 


-0.00030(0.00021) 



TABLE V: Simulation parameters and results obtained at L/a=16. 



f3 k Trajs. plq. 5t Acc. g| F M 

4.5500 0.1576500 32,309 0.540093(0.000014) 0.0143 0.785(0.003) 11.13131(1.41381) -0.00124(0.00018) 

TABLE VI: Simulation parameters and results obtained at L/a=18. 
C. improving discretization errors 

Since 0(a) discretization errors are not improved at all in the lattice actions, it is impor- 
tant to remove the scaling violation as much as possible. To do this, we perform the following 
improvements on the step scaling function and the DBF before taking the continuum limit. 

First let a(u,s) be the continuum limit of T, (u,s,l), i.e. a(u,s) = g§ F (sL) with u = 



1(3 



L/a 


N 


X 2 /dof 


a L/a,l 


a Li 


d.2 


a L/a,3 


a L/a,A a L/a,5 


6 


3 


9.0(1.3) 


0.4906( 


0.0025) 


-0.2749( 


0.0105) 


-0.1897( 0.0151) 




(3 


4 


1.4(0.5) 


0.5048( 


0.0014) 


-0.3993( 


0.0119) 


0.1136( 0.0283) 


-0.2042( 0.0184) 


6 


5 


1.3(0.5) 


0.5015( 


0.0032) 


-0.4240( 


0.0256) 


0.2538( 0.1301) 


-0.4043( 0.1815) 0.0899( 0.0808) 


8 


3 


2.1(0.7) 


0.5068( 


0.0018) 


-0.2308( 


0.0104) 


-0.2412( 0.0150) 




8 


4 


0.6(0.3) 


0.5153( 


0.0019) 


-0.3410( 


0.0260) 


0.0405( 0.0629) 


-0.1852( 0.0390) 


8 


5 


0.6(0.8) 


0.5153( 


0.0051) 


-0.3419( 


0.1697) 


0.0444( 0.7500) 


-0.1904( 0.9672) 0.0021( 0.3906) 


12 


3 


3.0(1.0) 


0.5239( 


0.0047) 


-0.1923( 


0.0118) 


-0.3019( 0.0198) 




12 


4 


1.1(0.6) 


0.5400( 


0.0038) 


-0.3614( 


0.0376) 


0.1063( 0.0884) 


-0.2671( 0.0550) 


12 


5 


1.0(0.6) 


0.5438( 


0.0039) 


-0.2783( 


0.0726) 


-0.2779( 0.2977) 0.2457( 0.3815) -0.2165( 0.1582) 


16 


3 


4.9(1.4) 


0.5308( 


0.0055) 


-0.1881( 


0.0266) 


-0.3057( 0.0375) 




1(3 


4 


1.8(0.8) 


0.5520( 


0.0039) 


-0.4948( 


0.0663) 


0.4387( 0.1516) 


-0.4762( 0.0903) 


16 


5 


1.9(0.9) 


0.5538( 


0.0050) 


-0.4403( 


0.1324) 


0.1801( 0.5648) 


-0.1283( 0.7332) -0.1457( 0.3025) 



TABLE VII: The results for the coefficients in the fit function ([2 

ggp(L). Its perturbative expression is given by 

a(u, s) = u + s Q u 2 + siu 3 + s 2 u 4 + ■ ■ • , (31) 

s = 6iln(s), (32) 

s 1 = \n(s)(b 1 2 \n(s) + b 2 ), (33) 

s 2 = ln(s) (V ln 2 (s) + h x b 2 ln(s) + b 3 ) , (34) 

where b^s are the coefficients of the /^-function introduced in sec. IIHI Recalling the para- 
metric form of the discretization error 



35| . the error normalized by a(u,s), denoted by 



So(u, s, I), is written as 



a(u, s) 



With eq. (1311) . the discretization error at the lowest order in u is found to be 

5 w (s,l) = (pi^-feiMs-O) - (pu-Mn(/)) =pi,..i-px,i-b 1 ]n(8)- ( 36 ) 

Now by replacing E (w,s,/) in eq. (1351) with H,x(u,s,l) = E (w,s,/)/(l + S^(s, I) u), the 
discretization error reduces to 0{u 2 ). Using Ei(it, s, I), the one-loop improved DBF is defined 

17 



by 



B 1 (u,s,l) = — - --. (37) 
Si(m, s, t) u 



This completes the one- loop improvement. 

The above procedure can be repeated to an arbitrarily higher order in u, but it requires 
the perturbative coefficients like pij and the perturbative expression of a(u, s) to the cor- 
responding order in u. All the coefficients necessary for the two-loop improvement are not 



available at this moment. Instead, we follow an alternative prescription proposed in Ref. I 
After the one- loop improvement, the scaling violation is written as 

= ExK s,l)- v(u,s) = s(2) u2 + Q 3 
cr(w, s) 

If one can somehow know <5^ 2 - ) (s, /), the scaling violation can be reduced to 0(u 3 ) by replacing 
E (w, s, Z) in eq. ( 1351) with 

E 2 (m, s, = E (u, s, /)/(l + 6 {1) (s, l)u + 5 (2) (s, I) u 2 ). (39) 

5^ 2 ^(s,/) can be determined by fitting our data for 6i(u,s,l) in eq. (1381) to the function 
quadratic in u. Notice that in order for this fitting to make sense, the perturbative series 
of a(u,s) must be known through 0(u 3 ). Since the first two coefficients, b\ and b 2 , are 
available, the correct value of a(u, s) can be calculated to 0(u 3 ) as seen from eq. (13T|) . 

5i(u, s, I) is fitted to the form of eq. fl38l) . neglecting 0(u 3 ) or higher order terms, for all 
possible pairs of (s, I) as shown in Fig. [2j The fit has to be performed in a weak coupling 
region where the perturbative expansion is reliable. We examine two fit ranges, < u < 1.6 
and < u < 2.0 to see the fit range dependence. The extracted values for 5^ 2 \s,l) are 
tabulated in Tab. IVIIII together with 6^(s, I) defined in eq. (J36J). 

The table shows that the values of 5^(s,l) and S^(s,l) lie between 10 -2 and 10 -3 , and 
S^ 2 \s,l) turns out not to depend on the fit range. In the following analysis, we employ 
8^ 2 \s,l) from the shorter fit range. It is also seen from the table that generally the coeffi- 
cients for (s,l) = (4/3, 12) are the smallest among others. This is anticipated because the 
improvement coefficient vanish as s approaches to unity or I becomes large. An exception is 
the one-loop coefficient 5 ( - 1 ^(4/3, 6). Since two-loop coefficient 5^ 2 - ) (4/3, 6) is, however, much 
larger than 5^(4/3,6), the smallness of 5^^(4/3,6) is probably by accident. In the data 
sets we have, the data with (s, I) = (4/3, 6) is the coarsest one. As we will show in the 
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{8,1) 


(4/3,6) (2,6) (8/3,6) (3/2,8) (2,8) (4/3,12) 




-0.00102 -0.0101 -0.0182 -0.00905 -0.0172 -0.00817 


5W(s,l) [0, 1.60] 
X 2 /dof 


0.0075(12) 0.0108(15) 0.0123(26) 0.0035(14) 0.0054(24) 0.0026(23) 
2.2 3.0 0.1 1.0 0.7 1.7 


S( 2 \s,l) [0, 2.0] 
X 2 /dof 


0.0061(9) 0.0106(12) 0.0112(18) 0.0038(13) 0.0053(18) 0.0014(17) 
2.2 2.1 0.2 0.9 0.6 1.4 



TABLE VIII: Coefficients for perturbative correction, 5^\s,l) and 5^ 2 \s,l), for each pair of (s, I). 
The square brackets in the first column indicate the fit range in u. 

following subsections, this data turns out to suffer from non-linear scaling violation larger 
than the linear one in the strong coupling region. Thus, we omit this data point throughout 
the analysis. Using 5^(s,l) thus obtained, we define the two-loop improved step scaling 
function ^{u, s, I) in eq. (J39l) . and in turn the two- loop improved DBF 

Ba(u,s,t) = 1 --. (40) 
h 2 [u, s,i) u 



D. strategy 

The continuum limit is taken for a fixed rescaling factor s and a fixed input length scale 
L varying a lattice spacing a (= L/l). As described in the preceding subsections, an input 
length scale is fixed by choosing a particular value of input coupling u. However, for a given 
s the number of data sets with different a in this work is, at most, two; (s, I) = (2, 6) and 
(2,8) for s — 2. While it is still possible to employ these two sets of data to evaluate the 
continuum limit, the validity of the linear extrapolation can not be tested. Alternatively, 
we may supplement a data set with a desired s by interpolating or extrapolating data of 
fi'sF^O' m ^- However, a lack of guiding principles in the interpolation or extrapolation may 
cause a systematic uncertainty. In this work, we use the two available data sets, (s, I) = (2, 6) 
and (2,8) to evaluate the continuum limit by linear extrapolation, and the other data sets 
are used to monitor the validity of the linear extrapolation. For this purpose, we introduce 
a relation which approximately converts the DBF for s' into that for s, as follows. 
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FIG. 2: Fit of 5\ to a quadratic function of u. The solid and dashed curves show the fit results 
and the fit ranges. 

We start with a closer look at the discretization error. The discretization error of the 
lattice DBF, i.e. BJu, s, I) — B(u, s) (z=0, 1,2), can be expressed in terms of an asymptotic 



expansion in 1/Z 35| as 

E>i(u, s, I) — B(u, s) 



I si 



ei (n) + 0(r 2 ), 



(41) 
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where e.i(u) is an unknown coefficient of 0(a) error and is a function of u. We then define 
the rescaled lattice DBF by 

B' t (u,s,l,s') = ^ Bl (u,s',l). (42) 

In addition, using the continuum counterpart of eq. (142 p . we define 

5B(u, s, s') = B{u, s) - B(u, s% (43) 

which represents the difference between the true continuum DBF and the rescaled continuum 
DBF. Combining eqs. ( T4T|) . P2|) and ( 143]) together and introducing 



ln(s') \l s'l 
we arrive at 

B'^u, s, Z, s') = £(u, s) + £(s, I, s') diu) - 5B(u, s, s') + 0(r 2 ). (45) 

Therefore, if SB(u, s, s') and 0(l~ 2 ) (or higher order) discretization errors are negligible 
compared to the statistical error of B[(u, s, I, s'), the numerical data of B'^u, s, I, s') plotted 
against £ will line up on a single line, and even two unknown coefficients in eq. (145]) . B(u, s) 
and ei(u), for given u and s can be extracted from that behavior. Instead, if both or one of 
them is large, the data will not align. Thus, whether B'^u, s, /, s') plotted against £ aligns 
or not tests the validity of the linear extrapolation within the statistical uncertainty. 

We comment on the size of 5B(u, s, s'). Solving eq. (CQ) perturbatively, the continuum 
DBF is found to be 

B(u,s) = -ln( S )|^ + M 2 ln( S )iM 2 + « 3 ln( S ) Q& a 2 & 2 ln(a) + hh + h 2 2 

+0(u 4 ln 2 (s)), (46) 

and thus the perturbative expression of 5B(u, s, s') is 

SB(u, s, s') = u 2 hx{s) In (^) - iftifta + u{ - ln(sa') - + ^ 2 2 ) 

+0 (w 4 ln(s) ln(s/s')) • (47) 



4 ), 



Since the numerical values of fej's are small, e.g. &i ~ 0.055, b 2 ~ —0.002, 6| F ~ O(10 
5B(u, s, s') is also small in the perturbative regime as 10~ 5 x u 2 (1.5 + 0.6u), 10~ 5 x u 2 (1.1 + 
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0.4m) and 10~ 5 x u 2 (-1.1 - 0.5 u) for (s, s') = (2, 16/12), (2, 12/8) and (s, s') = (2, 16/6), 
respectively. As u becomes large, 5B(u, s, s') may become sizable and at some point exceed 
the statistical error of B[{u, s,l, s'). Then, the alignment will be deformed. Notice that, 
the smaller B(u, s') is, the smaller 5B(u, s, s') is, and in particular, when B(u, s') = for a 
certain s', 5B(u, s, s') = holds exactly. 

We extract the continuum DBF B(u,s) as follows. First, we assume linear scaling vio- 
lation and calculate B(u,s) for s=2 by extrapolating the two data sets, (s',l) = (2,6) and 
(2,8), to £ = 0. Since s' = s, B[{u, s, I, s') = Bi(u, s', I) and 5B(u, s, s') = by construction. 
Thus we do not have to rely on the smallness of SB(u, s, s'). Then, to test the linearity of 
the scaling violation, we calculate the rescaled lattice DBF B'^u, s, I, s') with s=2 from the 
other data sets and plot them as a function of £(s, /, s'). If the data align within the statis- 
tical error of B'^u, 2, /, s'), the assumption of the linear scaling violation is valid, 5B(u, s, s') 
is negligible and then the value of B(u, s) thus obtained is reliable. Alternatively, once 
the linearity is confirmed, we can even determine the continuum limit by taking the linear 
extrapolation of B'^u, 2, 1, s'). Since SB(u, s, s') is negligible in perturbative regime, the lin- 
earity can be tested more rigorously in such a regime. When the data do not align, either 
or both of the linear violation dominance and small 8B(u, s, s') are invalid and the result 
for B(u, s) becomes uncertain. 

E. extraction of the continuum DBF 

Extrapolation to the continuum limit described in the following is carried out for every 
jack-knife ensemble, and the statistical error in the continuum limit is estimated by the 
single elimination jack-knife method. 

We begin with analysis at relatively weak coupling. Figure [3] shows the continuum limit 
of B'^u, s,l, s') for s = 2 (i=l, 2) at the four representative values of 1/u corresponding 
to u — 1.0,2.0,10/3,5.0, where the data with s' = 2 are shown in filled symbols and the 
other in open symbols and two data of the one-loop improved lattice DBF (B[) with s' = 2 
(filled squares) are linearly extrapolated to £ = 0. The two-loop improvement described in 
sec. IV CI is equivalent to tuning the improvement coefficients such that the resulting DBF 
reproduces the perturbative DBF in the region < u < 1.6. Indeed, the constant fit of 
the two-loop improved DBF with s' = 2 (filled diamonds) gives the value consistent with 
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FIG. 3: Linear extrapolation of B\ (filled squares) and constant fit of B' 2 (filled diamonds) to the 
continuum limit. The extrapolation and fit use the data with s' = 2 (filled symbols). The data 
with s' 7^ 2 (open symbols) are also shown to see whether they align or not. The values of (s', I) 
of the data shown are (4/3, 6), (2, 6), (8/3, 6), (3/2, 8), (2, 8), and (4/3, 12) from right to left. 
The data points are slightly shifted in horizontal direction for clarity. The perturbative predictions 
including the 2-loop (plus) and 3-loop (star) effects are also shown. 

the perturbative prediction when 1/u = 1.0 and 0.5 as seen in the figure. In the same 
region (l/w>0.5), the data of the one-loop improved DBF align and the linear extrapolation 
reproduces the perturbative DBF as well. Importantly, the extracted continuum DBF is 
clearly negative in this region. 

The deviation from the perturbative prediction appears at 1/u = 0.3, where the linear 
extrapolation gives the value closer to and consistent with zero. It is important to note that 
the data of the one-loop improved DBF align down to 1/u = 0.2 with a slope increasing with 
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u. From this observation, we conclude that, in the region 1/u > 0.2 (u < 5), SB(u, s, s') is 
small, the scaling violation is linear for the one-loop improved DBF and hence the extracted 
continuum limit is reliable. 

Next let us move on to the result at a stronger coupling shown in Fig. |H As seen from 
the figure, first the data of B[ except for the one with (s',/) = (4/3, 6) (right-most point) 
remains to align within the statistical uncertainty. Thus, the linear extrapolation of B\ is 
reliable at 1/u = 0.15. Secondly, the linear extrapolation of B\ and the constant fit of B 2 
lead to different continuum DBF. It appears that the constant fit of B2 is no longer valid 
and the linear fit appears to be more reasonable. Indeed, the linear fit of B 2 (solid line and 
open diamond at £ = 0) turns out to give the consistent limit as shown in the figure. 

From the alignment of B[, we infer that both SB and non-linear scaling violation re- 
main small. This is consistent with the fact that the continuum DBF obtained by linear 
extrapolation of B\ is consistent with zero and thus SB should be small as well. 

The deviation of the coarsest data from the linear behavior indicates that the linear 
discretization error no longer dominates others in the data with (s', /) = (4/3, 6). Since in 
general non-linear scaling violation can be large for small /, the data with I = 6 may suffer 
from this though it is not visible in the figure. To evaluate the potential uncertainty due 
to the 0(l~ 2 ) discretization error, we performed a linear fit without the / = 6 data. The fit 
result is shown as open square at £ = and the dashed line in Fig. HJ The result is consistent 
with that using the s' — 2 data only. 

From Figs. [3] and HI it turns out that for 1/u < 0-3 the extracted continuum DBF 
is consistent with zero. This indicates that in this region the running coupling constant 
reaches an infrared fixed point or, at least, the running appreciably slows down. In order 
to further investigate the existence of the infrared fixed point, we include the data obtained 
from 1=18 lattice at /3 = 4.55 into analysis. This data is combined with the data with / = 12 
to construct B[ with (s', /) = (3/2, 12). At = 4.55, the inverse SF coupling for 1=12 turns 
out to be 1/u = 0.107. On l=Q lattice, this value of 1/u is realized at /3 ~ 4.4. In such a 
small /3, the SF couplings are not calculated on 1=12, 16 lattices, and hence the following 
analysis is carried out without the data from 1=6 lattices. 

B[ constructed from the /=18 data is shown in Fig. [5] (filled circle). Since the four data 
points shown align well, we take the linear extrapolation using all of them and obtain the 
positive value in £ = 0. Interpretation of this result needs care as mentioned sec. [TT1 The most 
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FIG. 4: Same as Fig. [3l but the result at the stronger coupling is plotted. Solid lines denote 
the linear extrapolation using the data with s' = 2 (filled symbols). Dashed line shows the linear 
extrapolation using the data with the three smallest £. 

s=2.0, l/u=0.107 




FIG. 5: Same as Fig. U]but the data point obtained with (s', l)=(3/2, 12) (filled circle) is included 
in the analysis at l/u=0.107. The dashed line and the open square at £ = are the result of the 
linear fit. 



plausible explanation for this observation is that an IRFP exists in mirfp < 1.0/0.107 = 9.35. 



Figure [6] shows the 1/u dependence of the continuum DBF, where the results are compared 
with the perturbative calculations. It is seen that the running starts to slow down at 
around 1/u ~ 0.5, and eventually the coupling constant reaches a fixed point in the range 
of 0.107 < 1/u 0.3. When the DBF is positive, it is non-trivial for the continuum limit 
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FIG. 6: 1/u dependence of B(u,s) with s=2 obtained from the linear extrapolation of the data 
with s' = 2. Two- and three-loop perturbative predictions are shown by dashed and solid line, 
respectively. 

to exist. Thus we omit the positive DBF data from the figure. 




VI. SUMMARY AND OUTLOOK 



In this work, the running coupling constant of ten-flavor QCD is numerically investigated 
using lattice technique. The extrapolation of the DBF to the continuum limit is taken 
linearly assuming that the 0(a) scaling violation dominates the higher order ones. The 
DBF extrapolated approaches zero from below as the SF coupling constant u increases 
and when u>10/3 the DBF becomes consistent with zero. Further investigation at one 
particular strong coupling u = 9.3 (1/u = 0.107) is made using the data from the large 
lattice (I = 18), and suggests that the continuum DBF at this coupling is not negative. 
This indicates the existence of the infrared fixed point 10/3 < (?ir,fp ~ 9-3- The linear 
extrapolation is reasonably justified within the statistical error, but further rigorous check is 



clearly preferable. Combining our result with that of Ref. [12[, the critical number of flavors 
which separates the conformal phase and the broken phase is 8 < iVj rit < 10. 

In order to confirm the existence of IRFP or even determine the value of the fixed point 
more precisely, data from larger lattices with high statistics are necessary. It is, however, 
difficult to do with machines currently available to us, and probably more efficient methods 
or different approaches are necessary to go further. As mentioned in sec. [H the conformal 
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window can also be studied by looking at hadrons' spectroscopy or renormalization group 
analysis on the lattice. Currently the conclusions based on various methods are not consis- 
tent among them. In order to pin down N^ Tlt , these contradictions must be clarified with 
further studies. 

What is really important in the context of the WTC is the anomalous dimension of the 
ipip operator. The calculation of the anomalous dimension in ten-flavor QCD is on-going. 
The result will be published elsewhere. 

Once one has fixed an attractive candidate for WTC, the next important step would be the 



calculation of the S-parameter. The calculational method has been established in Ref. 



39|, 



where the QCD S'-parameter is calculated on the lattice for the first time and is correctly 



reproduced. Later, the method was applied to three-flavor QCD |40J, sextet QCD |4l| and 
six-flavor QCD 42J. In Ref. 42J, the evidence of the reduction of S'-parameter is reported. 



Another important quantity which should be calculated is obviously the mass spectrum of 
the candidate theory, including vector and scalar resonances, the decay constant of the NG 
boson and the chiral condensate. Although the precise determinations of these quantities are 
challenging, the direct comparison with the upcoming LHC results is extremely interesting 
and hence we believe that such calculations are worth a lot of efforts. 
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